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Abstract 

We examine the one dimensional Dirac equation with modulated or position dependent velocity. 
In particular, it is shown that using suitable velocity profiles it is possible to create bound state in 
continuum (BIC) like, as well as, discrete energy bound state solutions. 
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I. INTRODUCTION 



Since the rise of graphene l|, |2[ there have been a lot of studies on the application of Dirac 
equation in condensed matter systems. In such systems the effective low energy model of the 
quasi particles is described by the Dirac equation with a Fermi velocity vp which is 1/300 
times the velocity of light. The peculiar property of charge carriers in graphene is that they 
behave as relativistic massless fermions. Subsequently, there have been numerous studies on 
various theoretical as well as experimental aspects of graphene. In particular, confinement of 
the quasi particles using various types of magnetic and electrostatic fields have been studied 
{sj]. Other types of physically interesting states like the quasi bound states [4] have also been 
investigated. Here our interest lies in finding bound states in the continuum (BIC states) 
as well as discrete energy states in one dimensional heterostructures governed by the Dirac 
equation. The BIC are bound states embedded in the continuum [5|. Such states have 
a long history [6] and some years back the possibility of detecting BIC in semiconductor 
heterostructures has also been pointed out [7f. Later, BIC has been detected in photonic 
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been studied [10J]. 

On the other hand, there have been a number of studies on observable effects of modu- 
lated velocity or position dependent velocity in one dimensional heterostructures as well as 



graphene 
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14j . In such a scenario, the Fermi velocity is different in parts of the material. 



Backscattering in such systems described by one dimensional Dirac equation has recently 
been investigated thoroughly in ref Here our objective is to examine the(l + 1) di- 

mensional Dirac equation with position dependent velocity. To be more specific, we shall 
investigate the existence of BIC like states as well as discrete energy bound states in the 
presence of velocity profiles which continuously depends on the position. 



II. THE MODEL 



The Hamiltonian of a quasi-particle is of the form 

H = v F a x p x (1) 

where a x denotes the Pauli matrix. Now in the case of a heterostructure, it is quite possible 
that velocity may be different in two regions of the material. In other words, the velocity 
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becomes space dependent. However, this can not be achieved by simply replacing v F by v F (x) 
in Eq.fll]) since it would render the Hamiltonian non Hermitian. To keep the Hamiltonian 
Hermitian if the velocity vp becomes position dependent it is necessary to write it in the 
following way [ll] 

H = \Jv F (x)a x p x \Jv F {x) (2) 

It is easy to see that in case the velocity is a constant, the Hamiltonian (jSJ) reduces to the 
form p]). The Hamiltonian in Eq. (j2j) being hermitian we can also easily derive the associated 
continuity equation which expresses the conservation of probability: 

dp , dj x 



dt ~^ dx 



(3) 



where the current associated to the probability density p = ip^ip = ipt^Pi + ^2 ^2 is given by: 



jx = VvfiPI \/vFip2 + ex. 

The Hamiltonian in Eq. ()2]) operates on two component spinors ip(x) 
and the coupled equations for the components can now be written as 



y/v F (x) Vx V v f(x) ^{x) 



\/v F {x) p x \/vfJx) i>i{x) 



EMx) 



Eip 2 (x) 



(4) 

(ipi(x) , ip 2 (x)) T 

(5) 
(6) 



Now introducing an auxiliary spinor 111,2(2;) = ^/v F (x) ^1,2(2;), the above pair of equations 
may be written as 

v F p x ui = Eu 2 (7) 

v F p x u 2 = Eui (8) 

Let us now, for example, obtain the equation for the component U\ by eliminating u 2 from 
Eqs.© and © : 

, ri ,1,, 

; 2 = E 2 /h 2 (9) 



cPui 



+ v F v F 



e Mi, 



dx doc 

where the prime indicates differentiation w.r.t x. Let us now apply a further transformation 
involving the independent variable as 

1 



v F (t) 



dt 



and obtain from Eq.Q 



d 2 u\ o 



(10) 
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Eg. ( ITT]) looks like a free particle equation. However, whether or not it is a free particle 
equation depends on the range of z which in turn depends on Vp(x). We shall now examine 
the scenario with a couple of different choices of Vp(x). 



III. BIC LIKE SOLUTIONS 

Here we shall examine the existence of BIC like solutions of Eq. (fTTj) using different choices 
of the velocity vp(x). 
Example 1. 

Let us first consider 

vf(x) = vq cosh 2 (ax), v$ is a constant > (12) 

It follows that for a — > we recover the constant velocity scenario. Then from (fTOl) we 

find < z = tanh ax < and consequently Eq. ( fTTI) does not represent a free 
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particle. However, we choose a formal or unphysical solution [15[ of Eq. tfiTl) . namely, 

U\ = sin (ez) (13) 
Then the component ipi(x) of the original problem is given by 



il>i(x) = sech(aa;) sin I tanhaa; (14) 

\av J 

Clearly, the second part of ipi(x) is an oscillatory function bounded between ±1 while the 
first part is a function decreasing at both ends. Consequently, ipi(x) as given by Eq. (fl4|) is a 
normalizable solution for any value of the energy E and we conclude that Eq.( jT4"|) represents 
a BIC like solution. Now using the relation ([7]) the second component of the spinor can be 
found to be 



ip2(x) = —i sech(ax) cos 



tanh(ax) 



(15) 



av 

Clearly both components tpx^ are square integrable functions. One can indeed exactly eval- 
uate the normalization constant iV for the complete Dirac spinor. Choosing to normalize to 
1 the total probability density p = ip^ip we have: 

r+oo 

dxp(x) = / dx [\ipi(x)\ 2 + \Mx)\ 2 ] 



oc 



-N 2 (16) 
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FIG. 1: (Left panel) Plot of the component tpi (x), see Eq. [14] for the following set of values for the 
parameters: vo = 1, e = 10, a = 0.5 (solid curve), a = 1 (dashed curve), a = 2 (dot-dashed curve) 
(arbitrary units). (Right panel) Plot of the (normalized) probability density p{x) for the same 
values of the parameters as in the left panel. We clearly see that increasing the parameter a 
increases the localization of the BIC state. 

and the normalization constant is given by iV = y/a/2. We note that the two components of 
our spinor solution contribute differently, in general, to the total probability density. In any 
case, both components and consequently the Dirac spinor ip = (ipi,ip2) T represent a 

BIC like solution. In Fig. [T]we have drawn a plot of ip\(x) (left panel) and the normalized 
probability density p(x) = ift'if) = iplipi+ip^z (right panel) for a particular set of parameters 
in arbitrary units (see figure caption). It can be seen that the wavefunction oscillates with 
decreasing amplitude and eventually goes to zero. Also, the plot of the probability density 
(Fig. [T] right panel) shows that increasing the parameter a increases the localization of the 
BIC state. 

We can also easily check that the states described by Eqs. (114I15P have a vanishing 
probability current density as it is expected for a bound state: 

jx = Ny/v^ 



i) sin 



tanhax cos 



av 



tanh ax 



av 



(17) 
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Example 2. 

Here we choose 

Then it follows that 



2\2 



Vp{x) = v (l + ax ) 



a > 



7T 



2</av 



< z 



1 

2v~ 



x 



[1 + ax 2 ) y/a 



-L . — 

H — ■= arctan v ax 



< 



71 



2y/av 



Now proceeding as before the solutions iftipix) are found to be 



ipi{x) 



1 



1 + ax 2 
1 

i - 

1 + ax 2 



sm 



e 
2^ 



x 1 /— \ 

H — = arctan v«x 

1 + ax 2 y/a J 



cos 



e 

2^ 



1 + ax 2 



H — = arctan <Jax J 



(18) 



(19) 



(20) 



(21) 



As in the previous example both the components are square integrable functions and they 
contribute differently to the total probability density p = tp'ip. Normalizing to 1 the total 
probability density we find: 



dxp(x) 



dx [\Mx)\ 2 + \M%)\ 2 ] 



71 



N 2 



(22) 



where iV = a/ y/a/7T. 

It is readily seen that the solutions ( 1201) and ( J2TT) are valid for any value of e, and hence 
for any value of the energy eigenvalue E. In Fig. [2] we have plotted the component wave 
function ipi of Eq. (T2Q|) (left panel) and the normalized probability density (right panel) for 
the same set of parameters (in arbitrary units) as in Example 1. As in the first example, it 
is also seen here that the wave function oscillates with non constant amplitude before dying 
out. The probability density shown in the right panel again shows that with larger values 
of the parameter a there is stronger localization. 

From these two examples it is clear that one may choose many different profiles for v f (x) 
which would produce a BIC like solution. 



A. Bound state solutions 



In both the examples above we have considered unphysical solutions of Eq. (iTT|) . It now 
remains to be seen what happens when a physical solution i.e., a solution which vanishes 
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FIG. 2: (Left panel) Plot of the component tfji(x), see Eq. [20] relative to Example 2, for the following 
set of values for the parameters: vq = 1, e = 10, a = 0.5 (solid curve), a = 1 (dashed curve), a = 
2 (dot-dashed curve) (arbitrary units). (Right panel) Plot of the (normalized) probability density 
p{x) for the same values of the parameters as in the left panel. We clearly see also in this second 
example that increasing the parameter a increases the localization of the BIC state. 



at the boundaries is chosen. Here we shall only examine Example 1 in detail. In this case 
the problem is clearly that of a particle trapped in an infinite square well potential. The 
solutions of this problem are well known and are given by 



ui(z) 



sin 



irav n 1 
— = — (z + 



av 



a 2 v 2 n 2 n 2 



, n = 1,2,3, 



Hence the solutions of the original problem are given by 



ifti,n(x) = sech(ax) sin — (tanh ax + 1) 



7m 



havoim 



(23) 



(24) 



One may easily check that the solutions f l24|) are normalizable. Now again using relation (jjTJ) 
we obtain 

havo7m 



i)2,n( x ) — ~i sech(ax) cos — (tanh ax + 1) 

2 



Tin 



E- 



2.11 



(25) 



It may be observed that the solution ip2,n{x) represents the same energy level as ipi, n {%) 
Therefore the complete solution is given by 



sechax 




Tin 

— (tanh ax + 1) 

nn . 

— (tanh ax + 1) 

_ 



E„ 



havoTrn 



(26) 



So, in conclusion, whether or not the solutions of the original problem would represent a 
BIC or bound state depends on whether we choose an unphysical or a physical solution of 
the transformed equation given by Eq. liTTj) . 

IV. CONCLUSION 

Here we have examined the possibility of obtaining BIC as well as bound state solutions 
within the framework of one dimensional Dirac equation with a velocity which depends 
continuously on position. In particular, it has been shown that such solutions do exist for 
suitable velocity profiles. We feel it would be of interest to investigate BIC in the presence 
of a position dependent mass or a scalar potential as well looking for similar solutions in two 
dimensions. 

Acknowledgments 

One of us (P. R.) wishes to thank INFN Sezione di Perugia for supporting a visit during 
which this work was carried out. The Physics Department of the University of Perugia is 
also acknowledged for the kind hospitality. 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. 

Grigorieva, and A. A. Firsov, Science 306, 666 (2004). 
[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. 

V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005). 
[3] A. De Martino, L. Del? Anna, and R. Egger, Phys. Rev. Lett 98, 066802 (2007) 

A. De Martino, Solid State Comm. 144, 547 (2007) 

L. DellAnna and A. De Martino, Phys. Rev. B79, 045420 (2009). 

J. M. Pereira, F. M. Peeters and P. Vasilopoulos, Phys. Rev. B75, 125433 (2007) 

T. K. Ghosh, J. Phys: Condens. Matter 21, 045505 (2009). 

S. Kuru, J. M. Negro and L. M. Nieto, J. Phys: Condens.Matter 21, 455305 (2009). 
[4] PG. Silvestrov, Phys.Rev.Lett 98, 016802 (2007) 

A. Matulis and F.M. Peeters, Phys.Rev B77, 115423 (2008) 

8 



M. R. Masir, A. Matulis, and F. M. Peeters, Phys. Rev. B79, 155451 (2009) 
H. Chau Nguyen, M. Tien Hoang, and V. Lien Nguyen, Phys.Rev B79, 035411 (2009). 
[5] J. von Neumann and E. Wigner, Phys. Z. 30, 465 (1929). 

[6] J. Pappademos, U. Sukhatme and A. Pagnamenta, Phys.Rev A48, 3525 (1993) and references 
therein. 

[7] F. Capasso, C. Sirtori, J. Faist, D.L. Sivco, S.G. Chu and A.Y. Cho, Nature 358, 5656 (1992). 

[8] D.C. Marinica and A.G. Borisov, Phys. Rev. Lett. 100, 183902 (2008). 

[9] E. N. Bulgakov and A. F. Sadreev, Phys. Rev. B 78, 075105 (2008). 
[10] J.W. Gonzalez, M. Pacheko, L. Rosales and PA. Orellana, Euro.Phys.Lett 91, 66001 (2010). 
[11] N.M.R. Peres, J.Phys.Condens.Matter 21, 095501 (2009). 

[12] P.M. Krastajic and P Vasilopoulos, J.Phys.Condens.Matter 23, 135302 (2010). 

[13] A. Raoux, M. Polini, R. Asgari, A. R. Hamilton, R. Fazio, and A. H. MacDonald, Phys. Rev. 

B81, 073407 (2010). 
[14] A. Concha and Z. Tesanovic, Phys.Rev B82, 033413 (2010). 
[15] This solution does not satisfy the boundary condition at 



9 



